y[:,0] = X[:,0]**3+np.random.rand(n_samples)*20
linear_regression.fit(X, y) # train the model
y_predicted = linear_regression.predict(X)

plt.plot(X, y_predicted, 'k-')
plt.plot(X, y, 'yo')
plt.show()
